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In  a  recent  paper^,  the  problem  of  designing  a  feasible 
distillation  plant  for  the  production  of  heavy  water  was  dis¬ 
cussed.  Since  large  amounts  of  steam  are  required,  the  cost 
of  the  process  would  ordinarily  be  prohibitive  from  the  stand — 
point  of  plant  size,  and  the  quantity  of  fuel  required  for 

heating.  The  authors,  however,  have  in  mind  the  use  of  geo— 

2 

thermal  steam  ,  which  renders  the  essential  constraint  that  of 
plant  size. 

Under  various  assumptions  concerning  the  nature  of  the 
distillation  process,  the  problem  of  determining  the  most 
efficient  cascading  process  is  reduced  in  the  first  paper 
cited  above  to  a  multi-dimensional  maximization  problem,  which 
the  authors  solved  approximately  using  an  iterative  technique. 
Since  the  actual  physical  process  is  multi-stage,  it  may  be 
expected  that  the  theory  of  dynamic  programming  will  furnish 
a  more  systematic  computational  solution  to  questions  of  this 
type,  and  to  those  arising  from  more  realistic  assiwiptions .  In 
this  paper,  we  shall  consider  only  one  case  treated  by  the  cited 
authors . 
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Introduction 

In  a  recent  paper^,  the  problem  of  designing  a  feasible 
distillation  plant  for  the  production  of  heavy  water  was  dis¬ 
cussed.  Since  large  amounts  of  steam  are  required,  the  cost 
of  the  process  would  ordinarily  be  prohibitive  from  the  stand¬ 
point  of  plant  size,  and  the  quantity  of  fuel  required  for 

heating.  The  authors,  however,  have  in  mind  the  use  of  geo— 

2 

thermal  steam  ,  which  renders  the  essential  constraint  that  of 
plant  size. 

Under  various  assumptions  concerning  the  nature  of  the 

distillation  process,  the  problem  of  determining  the  most 

efficient  cascading  process  la  reduced  in  the  first  paper 

cited  above  to  a  raultl-dlmenslonal  maximization  problem,  which 

the  authors  solved  approximately  using  an  iterative  technique. 

Since  the  actual  physical  process  is  multi-stage,  it  may  be 
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expected  that  the  theory  of  dynamic  programming  will  furnish 
a  more  systematic  computational  solution  to  questions  of  this 
type,  and  to  those  arising  from  more  realistic  assumptions.  In 
this  paper,  we  shall  consider  only  one  case  treated  by  the  cited 
authors . 
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The  Analytic  Problem 


X.  Following  the  discussion  of  Karchetti^,  and  that  given  in 
A  X  mathematical  problem  is  that  of  minimizing  the^^plant  size ; 
for  an  ro-cascaae  process,  as  given  by 


J. 


subject  to  the  constraints 


m 

r  a. 


I j  a^  A  1  * 


Here  I  is  the  desired  total  enrichment,  and  a^  is  the 
enrichment  per  stage.  The  quantities  are  known  functions 
of  the  a^  determined  by  the  relations 


with  0  <  p  <  l5  Finally, 


I-P. 
^i  Pi 


N-j^  *•  1 . 5^10  i  Ng  =®  ^1^1^  ^i+1  *  ^i^i^  ^  *  ^t^t***i ®“1 , 


Vi  =  iNi- 


In  the  case  treated  in  ,  I  =*  500. 


^  C.  Marchetti,  Z.  Naturforschg .  9a,  1012  [195^3* 
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Setting 


(5) 


and  referring  to  (4),  wt  see  that  the  problem  is  equivalent  to 
that  of  minimizing  the  function 

V  »  g(ai )  +  — - —  -!■  ■  +  . . .  +  - r -  t  (6) 

1  a^  a^a^  ^1^2*-*^nt-l 

over  all  a^  subject  to  the  constraints 


(7) 


Variational  problems  of  this  type,  which  can  be  quite 
difficult  to  treat  by  conventional  methods,  can  be  resolved 
In  a  simple  fashion  computationally,  and  occasionally  analyti¬ 
cally,  using  the  techniques  of  dynamic  programming. 


Dynamic  Programming  Formulation 

Let  us  Introduce  the  sequence  of  functions,  defined 

as  follows: 

,,M  .  »»,  [.(.,)  .  ^  .  -jJ-  ,  - 

for  k  =  2,5,...,  where  the  a^  are  subject  to 


(9) 


and  X  may  assume  any  positive  value,  greater  than  one. 
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The  function  f^Cx)  is  readily  determined  from  (8)  and  (9), 


namely 


tAx)  «  Min  r  g(a^)  +  gCagVa^  1, 


(10) 


®1'®2 


where  =  x, 


Let  us  now  derive  a  recurrence  relation  connecting 


with  f.  (x).  Writing 


(*) 


.Min  [  g(a^)+  ^  [gCag)  +  S|22i 


+  •  •  »  + 


{"ij 


a^a^ 


we  see  that 


for  k  *  2,3, ..  . 


(12) 


This  is  an  application  of  the  principle  of  optimality, 
which  in  this  case  has  the  following  simple  physical  interpre¬ 
tation;  "whatever  enrichment  is  attained  in  the  first  cascade, 
at  whatever  cost  in  volume  of  plant,  the  remaining  enrichment 
is  to  be  obtained  using  minimum  plant  volume." 


The  solution  of  the  original  minimization  is  thus  reduced 
to  determining  the  sequence  j,  using  the  recurrence  re¬ 

lation  in  (12).  This  is  a  very  simple  process  which  can  be 
carried  out  via  a  direct  hand  computation,  a  direct  computation 
on  a  digital  computer,  or  by  using  variational  techniques. 


Discussion 
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The  technique  discussed  in  the  preceding  section  enables 
us  to  consider  more  realistic  processes.  If  we  allow  inhomo¬ 
geneous  cascades,  we  are  confronted  by  the  problem  of  minimizing 
an  expression  of  the  form 


ni  g.  {a. ) 

g,(a,)  +  Z  1  i 

i»2 


^i-l^i*“^m-l 


(13) 


where  the  sequence  ^g^(x)|is  loiown,  over  the  same  region  as 


above . 


To  treat  a  problem  of  this  type  introduce  the  sequence 
of  functions 


f  (x)  =  Min 

]  L 


a  •  T  a  •  • « .  a. 


i»k+l  “i-l“l‘ •  *“m-l 


(14) 


where  the  a^  are  subject  to 


Vk+r--^m  =  *' 


(15) 


for  k  =  1,2, .. . ,m— 1 . 


Then 


*  Min 

a  ,  ,a 
ra-1*  m 


®m— l^^m— 1^  ^  a 


®m^^m^  1 


itt-1 


(16) 


over  a  ^a  =«  x,  a  , ,  a„  >  1,  and  as  before, 
m— 1  m  m— i.  m  — 
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The  computational  solution  is  similar  to  that  for  the 
equation  of  (12). 


